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Abstract 

Considering the cascade decay D cC — > chB cbaA in which D,C,B,A are 
massive particles and c, b, a are massless particles, we determine for the shape of the 
distribution of the invariant mass of the three massless particles mate for the sub-set 
of decays in which the invariant mass ruab of the last two particles in the chain is (op- 
tionally) constrained to lie inside an arbitrary interval, mat G ['^afe* ^^'X]. 
An example of an experimentally important distribution of this kind is the "mg^/ 
threshold" - which is the distribution of the combined invariant mass of the visi- 
ble standard model particles radiated from the hypothesised decay of a squark to 
the lightest neutralino via successive two body decay,: q qx2 ~^ 0.^ qllXi^ 
which the experimenter requires additionally that mu be greater than mJ^"^/\/2. 
The location of the "foot" of this distribution is often used to constrain sparticle 
mass scales. The new results presented here permit the location of this foot to be 
better understood as the shape of the distribution is derived. The effects of varying 
the position of the mu cut(s) may now be seen more easily. 

Key words: Cascade decays, Kinematic Endpoints, Invariant Mass Distributions 
PACS: ll.80.Cr, 29.00.00, 45.50.-j 



1 Introduction 



We consider tlie cascade decay D cC cbB cbaA sliown in figure 1 
in which D, C, B, A are massive particles and c, b, a are treated as massless 
particles. 

In the context of the Large Hadron Collider (LHC), the most studied example 
of such a decay chain is the hypothesised decay of a squark to the lightest 
neutralino via three successive two body decays: q ~^ Q^l ~^ Q^^Xv The 

decay chain is relevant to more than just supersymmetric decays, however. 
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Fig. 1. The decay chain under study. Part (a) shows the chain in "generic" form 
with the massive particles on the "backbone" labelled with capital letters, and the 
massless visible particles labelled in lower case. Part (b) shows a commonly hypoth- 
esised squark decay chain open in many supersymmetric models. This particular 
chain is the example most frequently used to motivate the study of the "generic" 
chain in part (a). 

For example it may appear in models with extra dimensions, in which the 
particles would be Kaluza-Klein excitations. 

It has been frequently suggested [1,2,3,4,5] that one of the main ways we 
might hope to measure or constrain the masses of new particles produced at 
the LHC will be through the study of kinematic cndpoints of invariant masses 
constructed from the visible products of decay chains like those in figure 1. 

Given the observed four-momenta a^, and it is in principle possible to 
construct three non-trivial linearly-independent Lorentz-invariant quantities, 
namely m^^ = (a^ + b^)"^, ml^ = [h^ -\- c^Y ^ac = {^^ + ^^Y- ^ fourth 
quantity m^^^ = (a^ + h^^ + c^)^ is related to the first three via m^^^ = m^^ -|- 

In many models in which this decay chain is present, however, it is not pos- 
sible to tell two of the particles apart (typically a and h). For example, in 
the supersymmetric decay chain q — > ~^ ~^ Q^^Xi the particles a 
and b form a lepton-antilepton pair, and it may be impossible to say which 

observed Icpton was the "a" and which was the "&" . 1 In such models, the 
three independent variables are instead m^^, = max [m^^, m^^,] and 

(mg°™)^ = mm[m^^,ml^]. We see that the experimentally unidentifiable m^^ 
and have been replaced by experimentally observable {rrigi^^^)'^ and (m''^")^. 



If it is possible to measure the charge of the quark-jet some small fraction of 

the time, then there is a small discriminatory power in the a/b assignment to the 
observed leptons, if the spins of the sparticles are known. See [6] and subsequent 
papers. 
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Note that as before, m^^^ = m^^ + {rri^f^f + (mjj™)^. 

Although there are only three independent variables that categorise a given 
decay, by looking at more than one event of the given decay type it is often 
possible to construct more than three independent constraints on the masses 
of the particles involved. This has been done most frequently by constructing 
new variables which are arbitrary functions of the three independent variables, 
and then looking for the kinematic endpoints ^ of the distributions of these 
variables over a large number of events. 

In principle the shape of any one of these distributions contains more informa- 
tion than the position of the associated kinematic endpoint. Nonetheless, in 
previous years there has been a tendency to concentrate mainly on constraints 
coming from kinematic endpoints because interpreting these constraints does 
not require a detailed understanding of the way detector acceptance and recon- 
struction efficiency may vary across any of these distributions. Such detector 
effects may not be known with great confidence until a number of years after 
LHC turn-on. 

Recently, however, there has been a resurgence of interest in understanding 
the shapes of these kinematic distributions [7,8] not least because an under- 
standing of the shapes of these distributions is a pre-requisite for being able 
to fit the positions of the kinematic endpoints, even if one chooses to ignore 
all the information contained in the rest of the shape. In particular [8] has 
recently calculated the shapes of the distributions of the vriab, '^g/^^j ^"^^ 
rUabc distributions. One distribution whose shape was not calculated in [8] or 
elsewhere, however, was the shape of the m^^]^^^{x) distribution. 

Specifically the m^^^^^^'{x) distribution is a constrained form of the rUabc distri- 
bution in which the only events which are considered are those for which the 
value of ruab happens to exceed a cut x = nrf^ ™™ which is most frequently 
chosen to be 1/ \/2 times the maximum value which m^b can take. ^ The lower 
kinematic endpoint is the one whose position is most sensitive to the additional 
constraint (although the position of the upper endpoint can also be modified), 
and hence it is called a "threshold" rather than a "maximum". The ^^^^^^^{x) 
distribution is used because the position of the threshold, as a function of 
the masses of the particles involved, is more non-linear than the locations of 
other endpoints. This means that the position of the threshold can sometimes 
provide a handle on the absolute mass scale of the particles involved. Without 
a constraint from the endpoint of the vn}'^^^{x) distribution, the other edges 

^ i.e. the largest or smallest possible values that such a variable can take in a 
physically allowed decay configuration 

^ There is a long outstanding question (see remark in [3]) as to whether this some- 
what arbitrary choice for is indeed the one that gives the best constraint 
on the masses of the particles involved in the decays. 
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usually only succeed in constraining mass differences. 

Though important, the shape of threshold of the m*^^'^**'^(a;) distribution has 
been poorly studied in the past. Poor understanding of the shape of this 
edge has frequently resulted in the location of this endpoint being assigned a 
large systematic error (see for example [1,2,3,4,5]). Even allowing for a larger 
systematic error on the edge measurement, it often seems to be the case that 
the displacement of this edge from its "nominal" value is greater than a naive 
interpretation of the edge as a straight line would lead one to expect. 

We redress the historical imbalance in this note, by determining the shape of 
the m}^^^^{x) distribution in the spinless approximation, which is valid for su- 
persymmetry at the LHC. We do this for arbitrary values of a; = m^^* the 
imposed lower-bound requirement on niab- Given the form of the method, it 
turns out to be no extra effort to allow also an arbitrary upper-bound require- 
ment m'^^ '"^^^ on rriab, should one wish to exploit the information contained 
in the changes to the upper part of the distribution caused by varying this 
quantity. 



1.1 Existing results from earlier work 

First we note some existing results which we will make use of in later sections. 
[2] (and many of the references in the Supersymmetry Chapter therein) told 
us that: 




{ml - ml) {ml - m\) 
ml 



(1) 



Also [3,4] told us that 



'{ml - ml){ml - m\)/ml 
{ml — ml){ml — m\)/ml 
{mlml — mQm\){mQ — m^)/(m^m^) 
{mo - ruAf 



iff < mAfriD-i 
iff mAmjj < ml, 
iff mlmo < mAmQ, 
otherwise. 



(2) 



It was also shown in [3] and more recently in [9] that 




{x)r = x--i-2EP{Ei^^ 



(3) 



where 
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(p<c^,)^^^(<i±<^^ljlM^ and (5) 

c 

Ei?^^^ = V^^Tb^W (6) 
_ + (ml - m\) 

2mc ' 

where a; = rrf^ ™° is the smallest value of niab which events are allowed to 
have if they are to become part of the "threshold distribution" . 



2 Result 



The main result of this paper is that: 

1 dV 



(8) 



where 



Pirn 



abcJ 



f [J(/3) - J{a)] e{P-a) /Z if ml, < {mo - rriAf 
[ otherwise. 



In the above we have defined the function 

J{x) = — cosh~^ 

and defined the quantities 



ml + m\ — 
2mcmA 



(9) 



(10) 



/3 = min [Kf --^)^ (m+ f , (m-rf 

[Kf -^")^ K-,)^ 



Q; = max 



in which 



(11) 

(12) 



" ^ [^""^ + ^c)^L - (ml - ml) {{ml - m\) ^ Vx^}] (13) 



in which 



X = [mjj + - m„fej - AmArrij^ 



(14) 
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If the normalising constant Z in (9) is chosen to be 



Z = (^^^^^) («f - {jnt -^-f) , (15) 

then the integral of P{mlf,^) over physical values of m^f,^ is fixed at unity. This 
makes P(m^5^) into a normalised probability distribution, suitable for use in log 
likelihood fits, etc. The function ©(a;) in (9) is the Heaviside Step Function, 
equal to for a; < and equal to 1 for a; >= 0. If desired, the normalised 
distribution in mass space (rather than mass-squared space) is easily found 
via 

1 dV ^ 1 dr 
TTot d{mabc) Txot d{m^ahc) 



Remarks 



For the "classic" mg^rthreshold of [1,2,3,4,5], the following choices are made: 



«f '"^)' = WD', (17) 

(^cutmin)2^(^n.ax)2/2. (Ig) 



One can recover the standard "full" m^^,^ distribution (i.e. the distribution 
whose upper edge provides the measurement, not the constrained threshold 
distribution derived from it by applying extra cuts) by setting (mjj^* mm^2 
(18) to zero. 

Neither the final argument to the min in (11) nor the final argument to the 
max in (12) is strictly necessary These arguments arc there only to protect 
against accidental unphysical choices of (mj^^* max ^2 (mj^^* min-j2^ ^ g_ values 
which do not satisfy 

< (mS* '-'^f < (m^ "^n' < WD'- (19) 

The normalising constant Z in (15) has not been similarly protected, so care 
should be taken to ensure that (19) is satisfied when evaluating Z. 

If the function J{x) is only evaluated when (9) demands it, will always be 
non- negative and its square root may safely be taken in (13). 

Despite what the notation may appear to imply, either (or indeed both) of 
(m^j)^ and (^a^)' '^^^ t)e negative. It is important to retain the signs of these 
quantities when they are used in equations (11) and (12). 
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Fig. 2. Cartoon of the shape of the m^^^ distribution defined in (9). Note the dis- 
tribution is over the square of the mass, not over the mass itself. 

In mass-squared-space {i.e. plotting '^^ ^ rather than J the shape of 

^ abc' ^ abc) 

the distribution resembles a "top-hat distribution with sloping sides" as shown 
in figure 2. The flat upper part of the top-hat begins and ends at roots of 

^^cut max^2 _ (772+ )2 |-^cut mm^2 _ (772" )2 -^J^jf-Jj g^j^g ^J^g plaCCS at whlch 

the max and min functions in (11) and (12) have cusps. Though not straight 
lines, the "flanks" of the distribution are often not strongly curved, and so in 
the absence of effects from resolution, detector acceptance and other cuts, and 
in the absence of contamination from non-pure chains the lower edge should 
vanish linearly. This effect is not altered signiflcantly by the transition from 
mass-squared space to mass space, unless the threshold is itself very close to 
the origin, whereupon the edge takes on a quadratic character. 



3 Selection cuts and impurities 

Event selection cuts (and other effects) modify the shape of the mq//-threshold 
distribution. 

We illustrate this by generating a 100 inverse femtobarn sample of general 
SUSY events using Herwig 6.1 [10] and pass them through a simple detector 
simulation [11] which crudely models the response of the ATLAS LHC ex- 
periment. Following [1,2,3,4] and others, we work with the minimal SUGRA 
model known as "SUGRA Point 5" deflned by mo = 100 GeV, mi = 300 GeV, 
Ao = 300 GeV, tan/3 = 2.1 and /i > 0. This model has m^o = 121.5 GcV, 
mgfl = 157.2 GeV, m.^o = 233 GeV and has squarks of order 650 — 660 GeV 
(although the lightest stop is much lighter at 494 GeV). 

When constructing the "experimentally observed" mg^rthreshold distribution, 
we use exactly the procedure described in [3], which is to say in summary: 

• We require nieptons = 2, both leptons are opposite-sign same-family (OSSF) 
and > p?^ > 10 GeV. 

• We require njets > 2 and p^t —Pr — GeV. 

• We require -^t > 300 GeV. 

• We form only the larger of the two possible mgu combinations (one for 
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each of ji and J2) as our main interest is in preserving the quahty of 
the lower edge, which we do not wish to obscure with wrong-combination 
backgrounds. This jet-choice ensures wrong-choice backgrounds fall at high 
rather than at low mass values. 
• We repeat the above set of selection cuts but for opposite-sign different- 
family (OSDF) leptons. We add these events with weight -1 to the final 
histogram in order to subtract backgrounds from uncorrelated leptons. 

The resultant "experimentally observed" mg^-threshold distribution may be 
seen in figure 3(a). Figure 3(b) breaks this distribution down into the compo- 
nents in which the jet was correctly or incorrectly identified, and figure 3(c) 
shows the distribution predicted by the mass-space form of (9) using the 
masses m^o = 121.5, mj = 157.2, m^o = 233 and rrig = 656 GeV. The first 
three masses are taken directly from the model. The final mass represents the 
mass of a "typical" squark - an average of the masses of the twelve squarks 
participating in the chain. 

The figure shows that there is good agreement between the analytic shape 
in 3(c) and the "correct combination" component of the "experimentally ob- 
served" distribution, even though the jet-selection method should introduce a 
slight skew to the "correct combination" component. It seems that for these 
masses the skew is not very noticeable. The presence of the significant "wrong 
combination" background, however, means that it is important not to use the 
analytic shape on its own. 

When the m^H-threshold analysis was first proposed [1] the analytic form of 
the shape was not known. This motivated the deliberately high-mass-biasing 
jet-choice which seeks to keep the threshold of the mass distribution clean. 
Now that the shape of the distribution is known, it may prove more fruitful 
in later studies to bin both jet combinations, and then fit the analytic shape 
on top of a continuum background whose shape could perhaps be estimated 
from the side bands. 
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(a) (b) (c) 

Fig. 3. Three sets of mg;;-threshold distributions for the CMSSM mSUGRA model 
described in the text, (a) shows what might be seen in a typical LHC experiment 
after the cuts described in the text. Figure (b) shows (a) broken down into the four 
components which contribute to it. These arc (1) a narrow component coming from 
events in which the jet was correctly identified (solid line), (2) a broad component 
in which the jet was incorrectly identified (dashed line), (3) a small component 
representing the contribution from all SUSY sources which have nothing to do with 
the desired chain (dotted line), and finally (4) a small component from all SUSY 
sources which pass opposite-sign different-family versions of the cuts (dot-dashed 
line). The shape of component (4) matches the shape of (3) nicely and is added with 
negative weight in (a) to subtract component (3) . Figure (c) shows the shape of the 
m^//-threshold distribution as predicted in the text using the masses m^o = 121.5, 
mj = 157.2, m^o = 233 and ruq = 656 GeV. 
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4 Conclusions 



The shape of the distribution of the total invariant mass of the hght particles 
emitted in the decays shown in figure 1 has been calculated for the sub-set 
of decays which pass an additional upper and/or lower cut on the invariant 
mass of the last two particles emitted in the decay. This distribution is known 
generically as the mate-threshold distribution, or in its most commonly used 
specific example as the mg^-threshold distribution. Even though the calcula- 
tion treats all particles involved as scalars, the result is still exact in many 
cases, and will remain a good approximation in many others. The "classic" 
m^H-threshold distribution in supersymmetry is one of the places where the 
spinless approximation is exact, as spin effects from ql^l^ cancel those from 
ql~l^ . Deviations from the spinless approximation are potentially observable 
in UED and other non-susy models where spin effects in q*l~^l~ and q*l~l^ do 
not cancel completely, although the magnitude of these deviations are likely 
to be small as the production excess of q* over q* is itself expected to be small 
at the LHC over much of parameter space. 
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6 Appendix 

The derivation of the result (9) may be achieved using the following steps. 

• Work in the rest frame of the C. 

• In this frame the energy of the c is calculable (see (4)), as is the energy of 
the (a -|- 6)-system if it is taken to be an effective particle of known mass 
niab (see (6)). 

• The mass-squared m^^ may easily be shown to be distributed uniformly, 
and the polar angle 9 between the c and the (a + 6)-system in this frame 
will be uniform in cos 9 if either the spin of the C may be neglected or if we 
are unable (or choose not) to distinguish contributions from different spins. 
This will happen at the LHC where no distinction will be drawn between 
e'^e~ events and e~e'^ events. As a result we have: 




(20) 
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Change variables from cos 9 to m^^^, resulting in 

1 



r oc 



(21) 



The boundary of the region / may be shown in (m^^, m^j^)-space to be the 
intersection of the interior of one half of a hyperbola with a vertical strip of 



constant width in m„(, (starting at 



cut rmn\2 
ab 



) and ending at 



m 



cut max^2 
ab 



n 



It is possible to write the integration limits for I as upper and lower bounds 
on m^^ (respectively P and a in (11) and (12)) for given m^^^, allowing 
the order of integration to be exchanged, and the integral over m^^ to be 
performed analytically: 



roc 



1 



/ 



" '^rncPab 



—cosh 



(C)dimlb)di^lbc) 



m. 



c 



III 



A 



x=l3 



2 1 

abcJ ■ 



(22) 
(23) 



2mcmA 

from which the form of P(m^(,p) can be read directly. 
• The value of Z needed to normalise P{Tnlf^^) to unit area is found by per- 
forming the remaining integral. 
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